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Abstract 

Data  obtained  from  a  spatial  and  temporal  distri- 
bution study  of  sewage  effluent,  traced  with  fluorescent 
dye,  discharged  from  the  Municipality  of  Metropolitan 
Seattle  (METRO)  West  Point  treatment  facility  on  Puget 
Sound  are  presented. 

Following  an  initial  dilution  by  a  factor  of  about 
1/100  at  the  discharge  line  diffuser  ports,  effluent 
diluted  further  by  only  1/10  was  found  approximately 
6000  meters  north  (ebbing  current)  and  3500  meters  south 
(flooding  current).   Significant  concentrations  diluted 
by  a  factor  of  about  1/3  after  leaving  the  diffuser  ports 
were  found  near  shore,  and  concentrations  diluted  by  a 
factor  of  1/2.3  were  found  on  the  V,rest  Point  beach  (3.7 
ppt  effluent) . 

From  contours  of  sewage  concentration  and  density, 
wide  spatial  and  temporal  variations  were  observed  in 
both  throughout  the  5  days  that  the  experiments  were  con- 
ducted ,'30  August  through  3  September  1974).   Transverse 
cuts  of  the  plumes  indicated  that  contours  of  sewage  ef- 
fluent generally  coincide  with  density  contours  indi- 
cating that  density  stratification  is  influential  in 
confining  the  distribution  of  sewage  in  the  Sound.   A 
longitudinal  cut,  however,  indicated  significant  pene- 
tration of  density  layers  by  the  effluent  near  the  ends 
of  the  effluent  plume.   Implications  of  these  phenomenon 


vii 
on  modeling  effluent  diffusion  in  the  Sound  are  discussed, 
A  general  "buildup  of  dye  in  the  Sound  over  the  period  of 
the  experiments  was  not  observed. 

Comparison  of  these  dye  studies  with  drogue  studies 
carried  out  simultaneously  show  good  agreement  quanti- 
tatively for  the  30  August  data  and  qualitatively  for 
data  obtained  on  other  days. 


I .   Introduction 

From  30  August  to  3  September,  1974,  the  Applied 
Physics  Laboratory  (APL)  of  the  University  of  Washington 
undertook  a  study  of  the  sewage  effluent  distribution  in 
Puget  Sound  discharged  by  the  Municipality  of  Metropoli- 
tan Seattle  (METRO)  processing  plant  located  on  West  Point, 
Seattle.   This  study  addresses  two  principal  questions: 
Where  does  the  sewage  go?  What  is  the  nature  of  the  mix- 
ing process  acting  on  the  sewage  as  it  is  transported 
through^  Puget  Sound?  By  observing  the  nature  of  the  ans- 
wers to  these  two  questions,  conclusions  regarding  the 
applicability  of  model  types  can  be  drawn.   Whether  a 
buildup  of  sewage  concentration  occurs  or  whether  sewage 
residence  times  at  detectable  concentrations  are  limited 
is  investigated  as  well. 

Data  were  obtained  ~oy   cycling  a  towed  sensor,  "fish," 
through  depths  where  the  sewage  was  detected.   Fluores- 
cent dye  was  added  to  the  sewage  effluent  as  it  was  dis- 
charged from  the  processing  plant.   This  dye  was  detected 
by  the  fish  'as  it  was  towed  through  the  effluent  in  the 
Sound.   On  board,  real-time  instrumentation  allowed  ac- 
curate determination  of  dye  concentration,  fish  depth, 
and  vessel  position  with  respect  to  the  sewage  effluent 
discharge  ports  (outfall).   The  effluent  was  traced  from 
the  outfall  to  the  areas  in  the  Sound  that  it  reached. 
The  fish  was  capable  of  determining  temperature  and 
conductivity.   Density  data  determined  from  the  tempera- 
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ture  and  conductivity  data  are  also  used  in  this  analysis. 

Dye  concentration,  density,  depth  and  vessel  loca- 
tion data  were  recorded  on  magnetic  tape  on  board  and  then 
processed  ashore  to  provide  raw  information  used  to  hand 
contour  the  data.   The  visualizations  permitted  "by  the 
contoured  data  lead  to  most  of  the  observations  made  in 
this  report. 

An  investigation  of  current  properties  using  drogue 
information  was  conducted  by  Evans-Hamilton,  Inc.  in  paral- 
lel with  the  dye  studies.    Comments  concerning  the  con- 
clusions reached  by  these  studies  are  made  as  a  result  of 
the  data  obtained  from  the  dye  studies. 

This  investigator  'ore'nared  several  "orosrrams  for  use 
with  the  hand  held  computer-calculator,  HP-65,  for  navi- 
gation and  data  work-up  that  are  presented  in  Appendix  K. 

p 

II.   Experimental  Procedure 

Rhodamine  B  dye  was  injected  into  the  sewage  efflu- 
ent just  upstream  of  the  post-chlorinators  at  a  constant 
rate  of  1.73  kg/hr  (dry  weight)  for*   the  first  four  days 
of  the  experiment  (30  August  -  2  September).   The  flow 
from  the  outfall  (located  about  1100  meters  west  of  Vest 
Point  at  a  depth  of  73  meters)  was  70.9  million  gallons/ 
day  (mgd)  over  this  period  yielding  an  average  concentra- 
tion  of  dye  in  t^e  pipe  of  1.5  x  10   g/cc.   On  30  August 
the  average  flow  was  82.8  mgd  and  the  dye  injection  rate 
was  1.85  kg/hr  yielding  a  concentration  in  the  outfall 
of  1.4  x  10"'  g/cc.   A  dye  injection  rate  six  times  that 
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for  the  first  four  days  was  planned  for  the  fifth  day, 
but  cancelled  due  to  towing  cable  failure. 

Dye  concentration  was  measured  with  an  instrument 
package ( "fish")  towed  behind  APL ' s  50-foot  research  boat 
and  cycled  over  a  depth  appropriate  for  detecting  the  dye. 
Prior  to  pumping  dye  into  the  system,  a  background  fluor- 
escence level  equivalent  to  5  x  10"   g/cc  dye  concentra- 
tion was  determined  in  Puget  Sound  waters.   Since  dye 
concentration  in  the  outfall  is  about  1.5  x  10   g/cc, 
this  means  that  a  minimum  effluent  concentration  of  0.33 
parts  per  thousand  (ppt)  could  be  observed.   This  cor- 
responds to  a  maximum  dilution  factor  of  1/3000  obser- 
vable.  Towing  speeds  varied  from  3  to  6  knots.   The  fish 
consists  of  temperature,  electrical  conductivity,  dye  con- 
centration, and  depth  (pressure)  sensors  housed  in  a 
streamlined  aluminum  tube.   A  pump  within  the  fish  forces 
water  through  the  sensors.   An  electronics  package  accepts 
sensor  output,  multiplexes  the  data  onto  a  single  conduc- 
tor within  the  towing  cable  and  transmits  the  information 
to  shipboard  equipment.   Figure  1  shows  the  general  con- 
figuration of  the  fish  and  internal  equipment.   A  hydrau- 
lic winch  was  used  to  tow  and  cycle  the  fish  vertically. 
This  gasoline  powered  unit  is  capable  of  hauling  in  at  a 
rate  of  1  meter/second.   Both  water  depth  and  fish  depth 
indications  were  displayed  at  the  operators  station 
enabling  the  operator  to  cycle  the  fish  near  the  bottom 
while  avoiding  a  collision  with  it. 
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The  boat's  position  was  determined  to  within  a  few 

meters  using  a  Motorola  Mini-Ranger  radar  navigation 
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system  ,  which  determines  range  between  eaoh  of  two  on- 
shore transponder  units  and  the  shipboard  master  trans- 
mitting unit.   Ranges  were  continuously  displayed  on  the 
face  of  the  master  console  aboard  the  boat.   An  HP-65  pro- 
grammable calculator  was  used  to  convert  range-range  in- 
formation to  X-Y  coordinates  based  on  the  sewage  outfall 
position.   This  real  time  navigational  information  was 
essential  to  assure  that  the  boat's  navigation  closely 
followed  desired  paths  with  respect  to  the  outfall. 

Ranges,  dye  concentration,  fish  and  water  depths, 
and  clock  time  were  all  available  aboard  during  the  con- 
duct of  the  experiment.   Range,  all  fish  data,  water  depth, 
time,  and  manually  entered  "event  number"  were  recorded 
once  per  second  on  magnetic  tape  that  is  compatible  with 
APL's  CHI  2130  computer. 

Prior  to  the  experiment,  a  series  of  observations 
was  made  of  the  UW  Oceanography  Department's  Puget  Sound 
Model  in  order  to  provide  a  rough  feeling  for  the  expec- 
ted behavior  of  the  dye.   A  wide  degree  of  variability  from 
tidal  cycle  to  tidal  cycle  suggested  that  a  preset  pattern 
for  sampling  would  not  be  workable.   Instead,  it  was 
decided  that  a  series  of  arcs  crossing  the  Sound  in  east- 
west  directions  starting  near  the  outfall  and  working 
away  in  the  direction  of  the  current  would  be  the  most 
productive  means  for  detecting  the  area  containing  the 
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effluent.   In  addition  to  following  this  general  stra- 
tegy, sweeps  "back  up  the  plune(s)  were  conducted;  and 
several  excursions  through  water  outside  the  "boundaries 
of  well  defined  plumes  were  made  to  detect  detached 
pools  that  observations  of  the  Puget  Sound  Model  suggested 
might  exist. 

These  studies  and  the  drogue  studies  were  coordina- 
ted in  that  the  depth  at  which  dye  was  detected  was  com- 
municated to  Evans-Hamilton  personnel  for  use  in  setting 
drogue  depth,  and  the  establishment  of  flow  following 
slack  water  was  determined  from  the  drogues  and  used  to 
time  the  commencement  of  the  dye  detection  sweeps.   Nor- 
mallv  about  "50  minutes  were  allowed  to  T^ass  a"f*ter  ■flow 
commenced  in  order  to  assure  a  well  established  current 
pattern  prior  to  sampling. 

III.   Dye  Concentration  and  Density  Contouring 

Data  obtained  aboard  on  magnetic  tapes  was  returned 
to  the  APL  and  processed  on  the  CHI  2130  computer  to  plot 
dye  corcentration  and  density  profiles  for  each  cycle  of 
the  fish.   Figure  2  is  a  sample  of  this  output  that  was 
used  in  plotting  dye  and  density  contours.   The  boat's 
position  at  the  beginning  and  end  of  each  profile  was 
determined  from  a  printout  of  the  data  in  tabular  form 
as  shown  in  Figure  3.   Using  this  information,  plots  can 
be  made  that  correspond  to  the  X  and  Y  positions  the 
ship  traversed.   The  tracks  made  during  the  dye  detection 
sweeps,  however,  closely  approximated  straight  lines. 
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This  allowed  simplifying  the  contoured  data  to  be  plotted 
with  respect  to  distance  along  the  approximate  straight 
line  track.   The  straight  line  tracks  are  shown  as  the 
first  figure  of  each  of  Appendices  A  through  F.   Each 
plot  that  follows  in  these  Appendices  shows  distance 
along  the  respective  tracks  (corresponding  to  the  ap- 
propriate event  numbers)  with  the  zero  point  correspond- 
ing to  X  =  0  or  Y  =  0  depending  upon  whether  the  slope 
■  of  the  track  lies  within  4-5°  of  the  X  or  the.Y  axis  res- 
pectively.  Where  density  contours  have  been  plotted, 
they  were  done  on  translucent  paper  so  that  they  may  be 
overlaid  on  the  corresponding  dye -contours.   This  method 
greatly  improved  the  visualization  of  the  contoured  data 
and  leads  to  some  of  the  conclusions  expressed  in  this 
report.   The  depth  dimension  used  in  each  case  is  10 
meters/inch  with  the  water  surface  at  the  top  line.   All 
contoured  data  are  presented  as  though  the  observer  were 
looking  north  for  the  east-west  cuts  across  the  plumes, 
and  as  though  the  observer  were  looking  west  for  the  north- 
south  cuts.  '  Each  data  point  on  both  the  dye  concentra- 
tion and  density  contours  is  indicated  by  a  symbol  that 
represents  the  observed  concentration  at  the  point  plotted. 
Figure  4  provides  a  listing  of  all  symbols  used. 

The  reader  v/ill  note  upon  a  cursory  inspection  of 
the  plotted  data  that  it  is  closely  spaced  along  the 
path  of  the  fish  compared  to  the  horizontal  distance  bet- 
ween successive  depth  .cycles.   During  the  preliminary 
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Figure  4 
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stages  of  preparing  the  data,  computer  programs  that  were 
capable  of  contouring  were  investigated  and  it  was  found 
that  none  were  capable  of  handling  data  that  was  so 
unevenly  distributed.   Therefore,  hand  contouring  became 
necessary.   Although,  in  most  cases,  contours  could  be 
drawn  in  only  one  way,  occasionally,  the  data  became 
ambiguous  and  a  judgement  regarding  the  most  appropriate 
way  to  draw  the  contour  had  to  be  made. 

The  contoured  data  are  presented  in  six. segments 
in  Appendices  A  through  F.   These  correspond  to  the  data 
obtained  on  the  sweeps  indicated  in  Table  1. 

■ 

Table  1 

Experiment  Descriptions 

Day   Date   Flood/   North/     Direction  of    Appendix 
Ebb  '    South  of   Sweep  (Toward 
Outfall    or  Away  from 
Outfall) 

Away  A 

Away*  B 

Away  C 

Toward  D 

Single  loop        E 
made  in  area 

Single  loop        F 
made  in  area 


*  The  last  event  of  this  segment  was  a  sweep  back  up  the 
plume  toward  the  outfall. 


The  time  each  new  event  was  begun  is  indicated  on 


1 

30 

Aug. 

Flood 

South 

2 

31 

Aug. 

Ebb 

North 

2 

31 

Aug. 

Flood 

South 

2 

31 

Aug; 

Flood 

South 

h 

2 

Sep. 

Ebb 

North 

5 

3 

Sep. 

Ebb 

North 
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each  track  segment  on  the  track  plots.   The  reader  will 
note  the  significant  time  spans  encompassed  in  the  pat- 
terns presented.   It  is  obvious  that  the  data  presented 
is  not  synoptic.   Over  the  period  encompassed  during 
single  events,  the  contours  probably  do  approximate  a 
synoptic-  picture  in  that  during  the  short  period  of  time 
(about  30  minutes  more  or  less)  the  pattern  existing  in 
the  Sound  is  unlikely  to  be  changed  very  much.   Further- 
more, changes  that  do  take  place  within  this  time  frame 
are  likely  to  conform  to  the  random  processes  which,  pro- 
duce a  given  pattern  that  exists  at  any  given  instant. 
The  same  statement  is  not  observed  to  apply  over  time 
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Examples  of  significant  change  in  the  plume  patterns 
are  found  in  the  data  during  lacer  sweeps  of  the  plume 
area  within  the  same  tidal  cycle. 

The  times  of  the  slack  waters  that  preceded  and  fol- 
lowed the  plotted  tracks,  as  determined  from  tidal  current 
tables,  are  also  noted  on  each  track  plot. 

IV.   Observations 

The  data,  as  presented  in  Appendices  A  through  F, 
are  intended  to  be  useful  to  the  reader  for  the  purposes 
of  whatever  perspective  might  be  appropriate  to  his  in- 
terests.  Pertaining  to  the  specific  objectives  of  the 
METRO  study,  however,  some  observations  are  offered. 
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Variability  of  Sewage  Concentration  and  Density 

The  immediate  and  obvious  observation  upon  inspec- 
tion of  the  contoured  data  is  the  spatial  variability  of 
both  dye  concentration  and  density  that  was  encountered. 
Of  all  the  data  obtained,  day  1  contours  are  the  most 
uniform,  but  even  these  show  significant  variations  in 
depth  and  along  the  track.   Some  of  these  events  show 
two  areas  of  high  dye  concentration  instead  of  a  single, 
uniformly  developed  peak  (see  events  7  and  10.2,  day  1) 
that  might  be  expected  of  an  "idealized"  plume.   Intru- 
sions of  clean  water  (concentration  of  dye  less  than 
background)  are  seen  to  divide  the  plume  in  depth  also 
(see  events  11,  12,  and  13,  day  1). 

The  data  obtained  on  days  2,  4,  and  5  all  show  pat- 
terns of  dye  concentration  more  complex  than  that  of  day 
1.   Four  or  more  discrete  filaments  are  discernible  in 
much  of  this  data,  providing  an  indication  of  the  complex 
mixing  processes  occuring  in  the  Sound.   It  bears  mention- 
ing, again,  that  caution  must  be  exercised  in  analyzing 
the  contoured  data  presented  since  a  significant  amount 
of  time  passes  from  the  start  to  the  completion  of  each 
event. 

Models  and  Implications 

The  foregoing  observations  bear  significant  impact 
on  the  types  of  models  that  may  be  devised  to  describe 
the  effluent  distribution  in  the  Sound.   Many  models 
assume  that  the  diffusion  is  radially  symmetric,  and, 
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further,  that  probability  distributions  are  Gaussian  or 
of  a  shape  similar  to  that  of  a  Gaussian  distribution. 
From  the  foregoing  observations,  it  is  evident  that  nei- 
ther of  these  conditions  hold  in  the  Sound.   The  day  1 
data  most  closely  approximates  these  conditions,  but 
even  this  data  implies  that  such  models  must  be  used  in 
a  very  limited  fashion.   It  has  been  noted,  however,  that 
although  the  distributions  observed  hardly  conform  to 
the  conditions  mentioned  above,  a  good  fit  was  obtained 
on  day  1  for  a  model  that  assumes  the  "centerline"  con- 
centration, to  be  approximately  an  inverse  function  of 
distance  from  the  outfal"1  (concentration  proportional  to 
l/ya   Where  l<Ca^^/?^    The  drta  fro"  d?v  ^-    aD^sp*'  +^ 
conform  to  such  a  model  as  well,  but  only  3  events  were 
obtained  across  the  plume (s),  and  the  data  is,  therefore, 
not  conclusive.   Certainly,  any  model  that  attempts  to 
span  a  time  scale  greater  than  one  tide  could  not  suc- 
cessfully incorporate  the  assumptions  of  a  smooth  dis- 
tribution and  radial  symmetry. 

A  further  assumption,  that  the  effluent  may  be  des- 
cribed in  two  horizontal  dimensions,  is  frequently  made 
in  models  used  to  describe  flows  such  as  under  study  here. 
In  order  for  such  an  assumption  to  be  valid,  sufficient 
stratification  must  be  present  to  confine  the  effluent 
flow  along  density  contours  only.   It  will  be  seen  in  the 
following  observations  that,  although  stratification  greatly 
affects  the  distribution  of  effluent,  this  assumption, 
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also,  is  not  valid  during  portions  of  the  tidal  cycle. 

Depth  Distribution  of  Sffluent 

From  the  overlays  of  density  contours  it  is  seen 
that  the  stratification  in  the  Sound  is  influential  in 
the  observed  distribution  of  sewage  as  indicated  by  the 
dye.   Little  evidence  of  dye  contours  crossing  density 
contours  exists  in  the  east-west  crossings  of  the  plume(s) 
on  day  1.   Subsequently,  on  day  2  and  4  data  for  instance, 
the  trend  for  dye  contours  to  follow  density  remains  evi- 
dent, although  frequent  crossings  do  occur.   This  may 
be  partially  due  to  the  method  in  which  the  data  was 
processed.   Following  th^  first  day's  experiment,  a  sig- 
nificant lag  in  density  detector (s)  response  occurred, 
possibly  due  to  a  pump  or  piping  malfunction  (leak)  that 
went  undetected.   This  necessitated  averaging  the  density 
data  between  successive  half  cycles  to  eliminate  the  time 
lag.   Therefore,  the  contours  presented  for  days  2  and  4 
represent  values  averaged  over  from  50  to  200  meters, 
depending  upon  the  distance  covered  by  the  boat  between 
cycles.  .The  smaller  scale  density  fluctuations  were 
thus  averaged  out  while  the  dye  detection  equipment  con- 
tinued to  function  properly,  allowing  the  smaller  scale 
dye  concentration  structure  to  be  observed. 

It  is  noteworthy  that,  in  areas  very  close  to  the 
outfall  and  at  low  current  velocities,  the  effluent  was 
able  to  significantly  affect  the  density  structure. 
Event  5  on  day  1,  event  12.1  on  day  2,  and  event  2  on 
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day  4  all  show  a  divergence  of  density  contours  surround- 
ing areas  of  high  dye  (sewage)  concentration.   These  areas 
may  be  considered  "active"  in  the  sense  that  the  deeper 
portions  of  the  high  sewage  concentrations  are  lighter 
than  their  surroundings  and  the  shallower  portions  are 
heavier.-  These  conditions  were  generally  observed  within 
about  400  meters  or  less  of  the  outfall.   During  periods 
near  the  outfall  when  the  current  velocity  is  relatively 
low,  significant  vertical  penetration  of  the  dye  appears 
to  take  place  from  these  active  processes  and  may  lead  to 
the  observations  described  in  the  following  paragraphs. 

Longitudinal  Distribution  of  Sewage 

Event  11  on  day  2  is  a  longitudinal  track  from  north 
to  south  near  the  end  of  an  ebbing  current.   The  density 
contours  overlaid  upon"  the  dye  contours  clearly  show  that 
dye  penetrates  density  contours.   The  plume,  shown  to  be 
separated  by  about  1400  meters  in  the  portion  generated 
during  the  latter  half  of  the  tidal  cycle,  displays  ends 
in  each  segment  that  reach  a  lesser  density  layer  at  the 
northern  end  than  is  the  case  with  the  remainder  of  the 
segment.   Because  the  density  contours  are  penetrated, 
models  that  attempt  a  two  dimensional  description  cannot 
be  directly  applied.   (It  has  already  been  shown  that  a 
simple  plume  description  is  not  appropriate  due  to  the 
extreme  variability  of  the  flow  both  spatially  and  tem- 
porally.)  The  cause  for  this  observed  penetration  of 
the  density  layers  cannot  be  stated  with  much  confi- 
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dence  "because  of  the  limited  amount  of  data  available 
longitudinally.   Because  shallow  portions  of  the  plume 
segments  lie  at  the  ends,  it  is  possible  that  the  greater 
penetration  occurs  during  periods  of  low  current  vel- 
ocity that  occurs  during  a  change  of  tide  and  due  to 
changing  current  conditions  in  the  middle  of  a  tidal 
cycle  (perhaps  the  disestablishment  of  an  inshore  eddy 
significantly  slows  the  current  flow  over  the  outfall 
for  a  period  of  time  followed  by  a  return  to. greater 
velocities  as  the  new  trajectory  becomes  established) . 
During  this  period,  lateral  mixing  processes  might  be  more 
severely  restricted  and  thereby  diffuse  less  of  the 
buoyant  momentum  of  the  effluent  as  it  leaves  the  mixing 
ports  of  the  outfall.   This  would  tend  to  cause  the  ef- 
fluent to  penetrate  to  a  shallower  density  area  than  if 
the  current  velocity  were  higher.   Further  investigation 
is  needed  to  better  clarify  the  processes  that  cause  this 
phenomenon. 

Surface  Penetration 

It  is  significant  to  note  tha*o  no  dye  was  detected 
on  the  surface  during  these  experiments.   In  general, 
only  infrequently  did  dye  penetrate  to  above  20  meters, 
reaching  only  about  4  or  5  meters  at  the  most.   A  signi- 
ficant exception  to  this  occurred  on  day  5,  from  2000 
to  3500  meters  north  of  the  outfall.   Significant  con- 
centrations were  detected  to  about  4  meters  of  the  surface, 
and  it  is  likely  from  observations  of  the  contoured  data 
of  events  0  and  2.2  of  day  5,  that  dye  may  have  reached 
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a  depth  of  1  or  2  meters.   Significantly,  wind  condi- 
tions that  prevailed  during  day  5  were  different  (though 
not  particularly  severe)  than  on  the  previous  days  that 
the  experiments  were  conducted.   Wind  records  from  the 
West  Point  Coast  Guard  station  show  a  reversal  of  the  wind 
from  NE  on  day  4  to  SSE  on  day  5  with  an  increase  in 
average  speed  from  7.5  to  10.7  knots.   Although  the  mag- 
nitude of  this  speed  increase  does  not  seem  large,  a 
significant  change  of  surface  wave  conditions  were  noted 
on  day  5.   Also,  a  lesser  density  stratification  was  noted. 
These  conditions  obviously  had  an  impact  on  the  dye  dis- 
tribution that  was  observed.   It  would  be  expected  that 
an  even  greater  distribution  in  depth  might  be  noted  in 
the  winter,  and  this,  from  preliminary  indications  thus 
far  in  the  winter  experiment  analysis,  was  the  case. 
Dilution  of  Sewage  Effluent  in  the  Sound 
The  maximum  concentration  of  dye  observed  through- 
out the  conduct  of  the  experiments  occurred  during  a 
flood  tide  on  day  2,  event  12.1.   A  value  of  1.26  x  10   g/cc 
was  observed  400  meters  south  and  200  meters  east  of  the 
outfall.   This  represents  a  dilution  factor  of  1/119 
since  the  mean  concentration  in  the  outfall  on  day  2 
was  1.5  x  10   g/cc.   Of  interest  is  that  a  concentra- 
tion of  1.6  x  10~   g/cc  was  observed  over  6000  meters 
north  of  the  outfall  on  event  10,  day  2.   This  represents 
a  dilution  factor  of  only  1/7.9  in  the  distance.'   Since 
the  subsequent  event,  11  on  day  2,  shows  two  segments, 
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the  impression  that  the  northerly  segment  belongs  to 
a  previous  tidal  cycle  was  gained,  but  found  not  likely. 
Direct  measurement  of  current  velocity  by  drogues  was 
not  accomplished  on  this  day's  experiment.   However, 
by  using  the  maximum  current  observed  for  all  of  the 
drogue  studies  (day  1),  of  38  cm/sec,  and  the  fact  that 
the  average  current  during  a  tidal  current  of  the  dur- 
ation observed  is  approximately  .64  times  the  maximum 
velocity,  a  current  velocity  average  much  in  excess  of 
24.2  cm/sec  over  the  tidal  cycle  involved  is  not  likely. 
In  six  hours  at  this  speed,  a  distance  of  5230  meters 
could  have  been  traveled  by  the  leading  edge  of  the 
plume  without  taking  into  consideration  mixing  that 
occurs.   It  is  not  unreasonable,  therefore,  that  it 
reached  6000  meters  in  a  single  tidal  ebb  current.   Pos- 
sibly, during  the  portion  of  the  ebb  that  occurred  at 
relatively  high  current  velocities,  the  plume  was  directed 
toward  the  eastern  part  of  the  Sound  (toward  Shilshole  Bay) 
and  as  the  ebb  began  and  also  near  (but  not  at)  the 
end  when  it  slowed,  it  changed  to  a  more  central  route 
up  the  center  of  the  Sound.   Such  a  pattern  (and  the  many 
others  that  represent  such  widely  variable  character- 
istics) further  serves  to  emphasize  the  difficulty  in 
modeling  flow  in  the  Sound. 

From  these  observations,  it  is  evident  that,  as 
pointed  out  in  the  preliminary  report  to  METRO,  most  of 
the  mixing  and  dilution  of  the  effluent  occurs  at  the 
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diffuser  ports  at  the  outfall.-3  High  concentrations 
were  observed  south  of  the  outfall  as  well,  although  at 
a  distance  only  half  that  observed  to  the  north.   On  day 
2,  event  20,  a  concentration  of  4.0  x  10"'   g/cc  was 
observed  nearly  3000  meters  to  the  south  and  very  near 
shore.   This  concentration  represents  a  dilution  factor 
of  only  1/3.8  to  the  highest  concentration  observed  near 
the  outfall  on  this  day. 

The  data  indicates  that  sewage  reaches  considerable 
distances  north  and  south  of  the  outfall  and  very  near 
shore  with  little  dilution  from  mixing  in  the  Sound, 
within  one  tidal  cycle. 

■DU.X-LQ.UU  Ul   Ccwa.^c 

Throughout  the  five  days  the  experiments  were  con- 
ducted, a  buildup  of  a  dye  background  was  not  observed. 
All  dye  observed  showed  distinct  boundaries,  whether 
part  of  the  main  plume (s)  or  drifting  patches.   The  con- 
centrations within  the  plume  areas  remained  relatively 
consistent  and  clean  (less  than  background  dye  concen- 
tration) water  can  be  found  bounding  and  intruding  within 
the  main  concentrations.   Several  excursions  were  made 
outside  the  main  plume (s)  and  dye  was  not  detected.   It 
is  difficult  to  positively  identify  patches  of  dye  as 
belonging  to  a  previous  tidal  cycle,  although  event  4 
on  day  2  is  a  likely  candidate  since  the  location  and 
concentration  of  this  dye  patch  seems  to  indicate  that 
it  is  not  part  of  the  main  plume (s).   It  appears  that 
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most  of  the  sewage  that  enters  the  Sound  may  be  mixed 

to  below  background  concentrations  within  two  tidal  cycles. 

Sewage  Near  Shore 

Concentrations  of  effluent  ranging  as  high  as  2.7 
ppt  were  observed  on  several  occasions  very  near  shore, 
both  north  and  south  of  the  outfall.   This  is  somewhat 
evident  during  the  day  1  experiment  in  events  7  and  11, 
but  most  vivid  on  day  2  during  events  20  and  21  taken 
during  southerly  sweeps.   Observations  during  northerly 
sweeps  include  day  4  event  3.1,  and  day  5,  events  0  and  2.2 
As  shore  is  neared,  water  depth  decreases  and  the  depth 
to  which  the  fish  is  cycled  is  correspondingly  reduced 
and  this  is  seen  in  the  plots  of  the  above  tracks.   VJhere 
the  fish  depth  decreases  in  these  plots  significant  con- 
centrations of  dye  were  detected,  and,  in  most  cases, 
the  concentrations  actually  increase,  possibly  due  to  a 
near  shore  eddy.   That  effluent  reaches  shore  is  confirmed 
by  bottle  samples  taken  on  both  the  north  and  south  beach- 
es of  West  Point.   Dye  concentrations  as  high  as  5.1 
x  10"   g/cc  (corresponding  to  an  effluent  concentration 
of  3.4  PPt)  were  observed  on  the  beach  north  of  West 
Point  and  as  high  as  5.5  x  10"   g/cc  (3.7  ppt  effluent) 
south.   It  is  evident  that  sewage  diluted  by  a  factor  of 
1/2.3  of  that  found  to  be  the  maximum  after  exiting  the 
diffuser  ports  finds  its  way  ashore. 

V.   Amount  of  Dye  in  the  System 

In  an  attempt  to  determine  whether  significant  por- 
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tions  of  the  effluent  discharged  into  the  Sound  were 
going  undetected,  a  mass  balance  of  day  1  data  for  events 
10.2  through  12  was  made.   Day  1  was  chosen  since  the  plume 
pattern  was  more  uniform  that  day  and  a  good  cross  sec- 
tion of  the  dye  in  the  plume  was  obtained  for  each  of  the 
events  used.   Dye  contoured  data  was  averaged  for  each 
5  m  by  25  m  area  transverse  to  the  plume.   These  values 
were  then  added  to  provide  a  figure  for  the  total  amount 
of  dye  per  longitudinal  centimeter  as  demonstrated  by 
the  following  relation: 

cw[s/«n]=  (i  CiCggsFaJ  zsx.o'M  n 

where:   Cj_   =  the  concentration  observed  in  an  indi- 

6   P 
vidual  1.25  x  10   cm  cross  sectional 

area  (q/qwP/1.25   x  10   cm2) 

cev  =  a  linear  concentration  along  the  plume 

axis  (g/cm)  for  a  given  event  (ev) 

6   ? 
n   =  the  number  of  1.25  x  10   cm  squares 

for  which  a  concentration  (Cj_)  was 

observed  during  the  event 

The  values  of  cev  for  the  three  'events  were  averaged, 
thereby  entertaining  the  approximation  that  current  pro- 
perties were  the  same  as  each  of  these  segments  of  water 


passed  the  outfall.   Using  the  average  value  of  cev  = 
3.08  x  10    (g/cm) ,  and  assuming  that  this  average  fig- 
ure persisted  for  one. hour  (t  =  3600  sec),  and  further 
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assuming  that  this  segment  of  dye  passed  over  the  out- 
fall at  an  average  of  one  hour  following  the  commence- 
ment of  the  flood  tide  (an  approximation  based  on  the  time 
of  placement  of  the  drogues  into  the  water  and  the  observed 
position  of  the  boat  with  respect  to  the  drogues  during 
this  portion  of  the  experiment)  at  a  velocity  half  of 
the  maximum  observed  from  the  drogue  data  (v  =  .5  x  38 
cm/sec),  the  quantity  of  dye,  q  g/hr,  entering  the  Sound 
is  determined  as  follows: 


q  =  cev  x  v.  x  t 

=  (3.08  x  10"2)(0.5  x  38) (3600) 
s  2100  g/hr  (or  2.1  kg/hr) 


The  actual  injection  rate  measured  at  the  dye  pumping 
equipment  located  at  the  West  Point  sewage  treatment 
facility  was  1.73  kg/hr.   Agreement  between  these  two 
figures  is  good  and  there  is  no  indication  that  dye  is 
entering  the  system  and  going  undetected  into  areas 
outside  the  observed  plume. 

VI .   Comparisons  with  Drogue  Studies . 

The  data  obtained  on  day  1  were  digitized  and  ana- 
lyzed to  provide  statistics  that  could  be  directly  com- 
pared with  those  obtained  by  C.  C.  Ebbesmeyer  of  Evans- 
Hamilton,  Inc.  Density  contours  for  events  6  through  14 
of  day  one  were  used  to  prepare  digitized  data  for  each 
5  m  by  25  m  square  cross-section  (1.25  x  10   cm").   These 
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data  are  presented  in  Appendix  G-.   Values  shown  are  in 
tenths  of  a  power  of  ten  that  is  added  to  the  power  -10.3 
of  ten.   For  example;  a  value  of  4.5  shown  on  the  grid  is 
equivalent  to  a  power  of  0.45  to  be  added  to  -10.3, 
yielding  a  power  of  ten  equal  to  -9.85.   The  concentra- 
tion  represented  is  therefore  10"~~*  D   in  this  case,  or 
0.14  x  10   over  the  1.25  x  10   cm  area.   The  abscissa 
dimensions  are  arbitrarily  selected  digits  each  repre- 
senting 25  meters  along  the  track.   A  program  was  pre- 
pared for  the  HP-65  computer-calculator  to  calculate 
centroids  in  both  dimensions,  and  another  prepared  to 
calculate  the  standard  deviations  .and  variances  of  the 
distributions.   These  programs  are  included  as  Appendix  H. 

The  data  determined  in  this  way  are  presented  as 
Table  2.   A  plot  of  a  plume  that  corresponds  to  Figure 
3a  of  the  Evans -Rami It on  final  drogue  report'  is  pre- 
sented as  Figure  5.   Standard  deviation  data  were  converted 
to  represent  those  obtained  by  C.  C.  Ebbesmeyer  along 
the  axis  shown  in  his  Figure  7a.   A  plot  of  this  data 
versus  distance  from  the  outfall  is  presented  as  Figure  6, 
and  can  be  directly  compared  with  Figure  8  of  the  Evans- 
Hamilton  report  that  is  enclosed"  as  Figure  7  (overlay 
Figure  6  on  7  for  comparison). 

The  day  1  dye  experiment  analyzed  in  this  way  gen- 
erally confirms  the  drogue  data.   Differences  that  do 
occur  are  probably  more  a  result  of  the  difference  in 
time  that  the  dye  data  was  obtained  with  respect  to  that 
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Table  2 

Day  1  Dye  Statistics 


Event    X     Y       x         y         x         y 
No.     (m)    (n)   (104  cm)   (10^  cm)   (10°  cm)   (10°  cm) 
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53.5 

1.99 

10.8 

.  3.96 

1.18 

7 

183 

57.2 

2.21 

12.7 

4.91 

1.62 

10(1) 
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2.49 

10.0 

6.21 

1.00 

10(2) 

359 
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2.65 

15.0 
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2.42 

14.7 

5.87 

2.17 

12 

294 
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2.40 

16.2 

5.78 

2.61 

13 

411 
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3.78 

14.9 

14.31 

2.21 

14 

315 

56.4 
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of  the  drogue  data  (dye  data  is  about  g  to  1  hour  later 
than  the  drogue  data  oh  day  1).   Also,  a  qualitative 
inspection  of  data  from  experiments  conducted  on  subse- 
quent days  confirms  the  general  pattern  observed  using 
the  drogues.   However,  from  the  dye  studies,  it  is  seen 
that  distinct  filaments  are  observed  along  with  a  wide 
degree  of  variability  of  concentration  in  general,  and 
these  observations  would  indicate  that  the  current  pro- 
perties as  observed  using  drogue  data  are  likely  to  be 
significantly  optimistic  with  respect  to  the  concentra- 
tions that  they  predict.   It  is  not  possible  to  detect 
individual  filaments  using  drogues,  and  it  is  evident 
that  mu  1  ~t i  *"f/3 1  3 men^*   TDiunies  t^qvct  g+  tlprou r~,^~' c^i ■? i *?■.  +.Vi p>  ^qtq 
Eddy  coefficients  therefore  cannot  generally  be  predicted 
confidently  using  models  consistent  with  the  smoothly 
averaged  type  of  data  obtainable  from  drogues. 

VII.   Summary 

Most  of  the  mixing  over  a  single  tidal  cycle  ap- 
parently takes  place  at  the  diffuser  port  at  the  outfall. 
The  maximum  concentration  found  in  the  vicinity  of  the 
outfall  was  1.26  x  10   g/cc,  representing  a  dilution  of 
the  sewage  by  a  factor  of  1/119.   After  leaving  the  vici- 
nity of  the  outfall,  very  large  spatial  and  temporal 
variations  of  sewage  concentrations  are  found  in  the 
Sound.   (The  same  can  be  said  of  density  profiles.)   In- 
dications of  sewage  diluted  by  a  factor  of  less  than 
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l/lO  in  the  Sound  proper  wore  found  as  far  as  6  km  north 
and  3.5  km  south  of  the  outfall.   Near  shore,  concentra- 
tions as  high  as  4  x  10    s/cc  were  noted,  representing 
a  dilution  by  a  factor  of  1/3.1  since  entering  the  Sound. 
Sample  "bottles  ashore  were  found  to  contain  concentra- 
tions as  high  as  5.5  x  10"   g/cc,  representing  a  dilution 
by  a  factor  of  just  1/2.3. 

The  data  indicate  that  the  density  stratification  is 
influential  in  confining  the  distribution  of  sewage  in 
depth,  and  little  significant  crossing  of  dye  and  den- 
sity contours  occurs  in  the  transverse  cuts  across  the 
plumes.   A  longitudinal  cut  back  up  a  plume  near  the  end 
of  an  ebb  current,  however,  demonstrated  significant  cross- 
ing of  density  contours  by  dye  concentration  contours 
which  was  most  pronounced  near  the  ends  of  the  plumes. 
This  may  indicate  that  the  sewage  is  more  effective  in 
penetrating  the  pycnocline  over  the  outfall  during  low 
current  velocities. 

Models  that  are  used  to  describe  diffusion  of  the 
nature  occuring  from  the  METRO  sewage  outfall  at  V,rest 
Point  often  assume  that  the  flow  is  two  dimensional,  or 
even  that  the  diffusion  is  radially  symmetric.   Experi- 
mental observations  indicate  that  none  of  these  assumptions 
can  be  applied  to  the  conditions  existing  in  the  Sound 
without  accepting  very  significant  error  except  over  an 
extremely  limited  range  of  space  and  time  scales.   If 
two-dimensional  models  are  employed,  they  should  be  ap- 
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plied  to  that  portion  of  the  tidal  cycles  that  contain 

well  established  currents,  since  some  evidence  exists 
that  sewage  penetrates  density  contours  more  during  the 
initial  stages  of  the  cycle.   A  model  that  assumes  the 
"centerline"  concentration  to  be  approximately  an  inverse 
function-  of  distance  from  the  outfall  was  applied  to 
the  data  with  limited  success.   Such  a  model  applies 
during  tidal  cycles  that  produce  a  plume  that  is  well 
defined  and  contains  little  transverse  variation.   Such 
a  condition  leads  to  less  mixing  than  in  plumes  contain- 
ing more  variability  (ie:  more  filaments  or  eddies). 
Only  the  data  from  day  1  provided  a  significant  indica- 
tion for  success  of  such  a  model. 

The  buildup  of  a  general  dye  background  was  not  ob- 
served although  it  is  not  possible  to  say  that  the  clear 
water  observed  in  the  experiment  was  a  result  of  advec- 
tion  of  dye  out  of  the  observation  range  (unlikely) ,  or 
a  result  of  dye  having  been  diluted  to  below  the  limits 
of  detectability.   This  experiment  did  not  address  the 
question  of  what  happens  after  a  tidal  cycle  is  complete 
with  respect  to  whether  dye  is  transported  back  across 
the  outfall,  thereby  surviving  the  current  reversal. 
The  winter  experiment,  under  analysis  at  this  time,  in- 
dicates that  significant  patches  of  dye  do,  in  fact, 
survive  the  reversal  in  current  and  are  transported  back 

o 

across  the  outfall. 

Comparison  of  the  dye  studies  data  with  drogue  stu- 
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dies  data  indicates  that  the  two  are  generally  consistent. 
Significantly,  however,  the  variability  of  conditions 
that  were  found  in  the  dye  studies  cannot  be  determined 
by  drogue  studies.   Therefore,  the  drogue  data  probably 
underestimates  the  peak  concentrations  likely  to  be 
observed -in  nature. 
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Appendix  A 


Day    Date    Flood/    North/      Direction  of 

Ebb  '     South  of    Sweep  (Toward 
Outfall     or  Away  from 
Outfall) 

1    30  Aug.   Flood     South         Away 


Event  No. 

Type  of  Contour 

Page  No. 

(Track  Chart 

,  .   35) 

05 

Density 

36 

05 

Dye 

37 

06 

Density 

38 

06 

Dye 

39 

07 

Density 

40 

07 

Dye 

41 

10(1) 

Density 

42 

10(1) 

Dye 

43 

10(2) 

Density 

44 

10(2) 

Dye 

45 

,-.1 

Density 

46 

11 

Dye 

.   47 

12 

Density 

48 

12 

Dye 

49 

13 

Density 

50 

13 

Dye 

51 

14 

Dye 

52 

53 

Appendix  B 


Day    Date    Flood/    North/      Direction  of 

Ebb  '     South  of    Sweep  (Toward 
Outfall     or  Away  from 
Out fell) 

2    31  Aug.   Ebb       North  Away* 


**  The  last  event  of  this  segment  was  a'  sweep  back  up 
the  plume  toward  the  outfall. 


Event  No. 

Type  of  Contour 

Page  Mo. 

(Track  Chart 

.   54) 
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Dye 
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Dye 
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Dye 
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Dye 
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Dye 
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Density 
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Dye 

61 
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Density 
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Dye 
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Dye 
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Appendix  D 


Day 


Date 


Flood/ 
Ebb 


North/ 
South  of 
Outfall 


Direction  of 
Sweep  (Toward 
or  Away  from 
Outfall) 


31  Aug.   Flood 


South 


Toward 


Event  No. 


Type  of  Contour 


Page  No. 


(Track  Chart  74) 

20  Dye  75 

21  Dye    •  76 

22  Dye  77 
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Appendix  E 


Day    Date    Flood/    North/      Direction  of 

Ebb      South  of    Sweep  (Toward 
Outfall     or  Away  from 
Outfall 

4     2  Sep.   Ebb       North       Single  loop 

made  in  area 


Event  No. Type  of  Contour Pa^e  Ho 

(Track  Chart  79) 

1  Dye  80 
2m                          Density           81 

2  Dyr  82 

^  ./ 

84 
85 


3.1 

Dye 

3.2 

Dye 

4.2 

Dye 
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Appendix  F 


Day    Date    Flood/    North/      Direction  of 

Ebb      South  of    Sweep  (Toward 
Outfall     or  Away  from 
Outfall 

5     3  Sep.   Ebb       North       Single  loop 

made  in  area 


Event  No. Type  of  Contour Pap;e  No . 

(Track  Chart  87) 

0  Dye  88 

2.1  Dye  89 

2.2  Dye  90 
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Appendix  G- 
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Appendix  H 
HP-65  Programs 
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